Plasmon excitations and ID - 2D dimensional crossover in quantum crossbars 
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Spectrum of boson fields and two-point correlators are analyzed in quantum crossbars (QCBs, a 
superlattice formed by m crossed interacting arrays of quantum wires), with short range inter- wire 
capacitive interaction. Spectral and correlation properties of double (m = 2) and triple (m — 3) 
QCBs are studied. It is shown that the standard bosonization procedure is valid, and the system 
behaves as a sliding Luttinger liquid in the infrared limit, but the high frequency spectral and 
correlation characteristics have either ID or 2D nature depending on the direction of the wave 
vector in the 2D elementary cell of reciprocal lattice. As a result, the crossover from ID to 2D 
regime may be experimentally observed. It manifests itself as appearance of additional peaks of 
optical absorption, non-zero transverse space correlators and periodic energy transfer between arrays 
("Rabi oscillations"). 
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I. INTRODUCTION 

The behavior of electrons in arrays of ID quantum 
wires was recognized a challenging problem soon after the 
consistent theory of elementary excitations and correla- 
tions in a Luttinger liquid (LL) of interacting electrons in 
one dimension was formulated (seeu for a review). One 
of fascinating targets is a search for LL features in higher 
dimensionscl. Although the Fermi liquid state seems to 
be rather robust for D > 1, the possible way to retain 
some ID excitation modes in 2D and even 3D systems 
is to consider highly anisotropic objects, in which the 
electron motion is spatially confined in major part of the 
real space (e.g., it is confined to separate linear regions 
by potential relief). One may hope that in this case weak 
enough interaction does not violate the generic long-wave 
properties of the LL state. Arrays of interacting quantum 
wires may be formed in organic materials and in striped 
phases of doped transition metal oxides. Artificially fab- 
ricated structures with controllable configurations of ar- 
rays and variable interactions are available now due to 
recent achievements in nanotechnology (see, e.g., Refs. 

ID- 

The simplest 2D ensemble of ID nanoobjects is an 
array of parallel quantum wires. The conventional LL 
regime in a single ID quantum wire is characterised by 
bosonic fields describing charge and spin modes. We con- 
fine our discussion to the charge sector (LL in the spin- 
gapped phase) . The Hamiltonian of an isolated quantum 
wire may then be represented in a canonical form 



L/2 
-L/2 



(1) 



Here L is the wire length, v is the Fermi velocity, 8, ir 
are the conventional canonically conjugated boson fields 
and g is the dimensionless parameter which describes_i-he 
strength of the interaction within the chain (see, e.g., Bel). 
The interwire interaction may transform the LL state ex- 
isting in isolated quantum wires into various phases of 



2D quantum liquid. The most drastic transformation is 
caused by the interwire tunneling t± in arrays of quantum 
wires with intrawire Coulomb repulsion. This coupling 
constant rescales towards higher values for strong interac- 
tion (g < 1/2), and the electrons in array transform into 
2D Fermi liquids. The reason for this instability is the 
orthogonality catastrophe, i.e. the infrared divergence in 
the low-energy excitation spectrum that accompanies the 
interwire hopping processes. 

Unlike interwire tunneling, the density-density or 
current-current interwire interactions do not modify the 
low-energy behavior of quantum arrays under n certain 
conditions. In particular, it was shown recentlyuEl that 
an interaction of the type W(n — n'), which depends on 
the distance between the wires n and n' but does not con- 
tain current coordinates x,x', imparts the properties of 
a sliding phase to 2D array of ID quantum wires. In this 
state an additional interwire coupling leaves the fixed- 
point action invariant under the " sliding" transformation 
n — ► n + a n and 7r n — * 7r„ + a' n . The contribution of 
interwire coupling reduces to a renormalization of the 
parameters v — > v(q±), g —* g{q\_) in the LL Hamilto- 
nian (]]]), where q± is a momentum perpendicular to the 
chain orientation. Such LL structure can be interpreted 
as a quantusj analog of classical sliding,phases of coupled 
XY chainsE2l. Recently, it was founcO that a hierarchy 
of quantum Hall states emerges in sliding phases when a 
quantizing magnetic field is applied to an array. 

In the present paper we concentrate on another aspect 
of the problem of interacting quantum wires. Instead of 
studying the conditions under which the LL behavior is 
preserved in spite of interwire interaction, we consider sit- 
uations where the dimensional crossover from ID to 2D 
occurs. Dimensional crossover is quite well studied e.g. 
in thin semiconducting or superconducting films where 
the film thickness is the control parameter that rules 
the crossover (see e.g. Ref. 12. It occurs in strongly 
anisotropic— systems like quasi-pne-dimensional organic 
conductors^ or layered metalalj. In the latter cases 
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temperature serves as a control parameter and crossover 
manifests itself in interlayer transport. In metals the lay- 
ers appear "isolated" at high temperature, but become 
connected at low temperatures to manifest 3D conduct- 
ing properties. Here we intend to study another type of 
dimensional crossover, i.e. a geometrical crossover, where 
the phase variable serves as a control parameter, and the 
excitations in quantum array demonstrate either ID or 
2D behavior in different parts of reciprocal space. 

The most promising type of artificial structures where 
this effect is expected is a periodic 2D system of m cross- 
ing arrays of parallel quantum wires. We call it "quantum 
crossbars" (QCB). The square grids of this type consist- 
ing of 2 arrays ware-considered in various physical context 
in early paperal3il3. In Refs. [l7],[l^ the fragility of the 
LL state against interwire tunneling in the crossing areas 
of QCB was studied. It was found that a new periodic- 
ity imposed by the interwire hopping term results in the 
appearance of a low-energy cutoff A; ~ hv/a where a 
is a period of the quantum grid. Below this energy, the 
system is "frozen" in its lowest one-electron state. As a 
result, the LL state remains robust against orthogonality 
catastrophe, and the Fermi surface conserves its ID char- 
acter in the corresponding parts of the 2D Brilllouin zone 
(BZ). This cutoff energy tends to zero at the points where 
the one-electron energies for two perpendicular arrays 
and 6k 2 become degenerate. As a result, a dimensional 
crossover from ID to 2D Fermi surface (or from LL to 
FL behavior) arises around the points ep t = ep 2 . 

We study this dimensional crossover for Bose excita- 
tions (plasmons) described by canonical variables 9, n in 
QCB. In order to unravel the pertinent physics we con- 
sider a grid with short-range capacitive inter-wire inter- 
action. This approximation seems natural for 2D grids 
of carbon nanotubestl, or artificially fabricated bars of 
quantum wires with grid periods which exceed the lat- 
tice spacing of a single wire or the diameter of a nan- 
otube. It will be shown below that this interaction can 
be made effectively weak. Therefore, QCB retains the 
ID LL character for motion along the wires similarly to 
the case considered in Rcf. |[ At the same time, the bo- 
son mode propagation along some resonant directions is 
also feasible. This is essentially a 2D process in the 2D 
BZ (or in the elementary cell of the reciprocal lattice). 



We start the studies of QCB with a do uble QCB m 
(section ||). In the first two subsections [I A and 



II B 



= 2 
we 

introduce basic notions and construct the Hamiltonian 
of the QC B. T he main approximations are discussed in 
subsection II C. Here we substantiate the used method 



(separable interaction approximation) and show that in- 
teraction between arrays in QCB is weak. The energy 
spectra for square QCB an d til ted Q CB are describ ed in 
detail in two parts II D l| and IIP 2] of subsection [ID. 
Various correlation functions and related experimentally 
observable quantities (optical absorption, space correla- 



tors) are discussed in the last subsection II E . We predict 
here effect of peculiar "Rabi oscillations" - periodic en- 
ergy transfer from one of the QCB array to another. 

Triple QCB (to = 3) formed by three arrays lying in 
parallel planes are studied in Section III. Such hexago- 
nal grida.may be useful for three-terminal nanoelectronic 
deviceso. The plasmon spectra of triple QCB possess 
some specific features in comparison with double QCB. 
We introduce the main notions and construct the Hamil- 



tonian of symmetric triple QCB (subsection III A ) , anal- 
yse the peculiarities of the frequency spectrum (subsec- 
and illustrate them by description of triple 
perio dic en ergy transfer between all 



[II B) 



tion 

Rabi oscillations 
three arrays (subsection QIC). The results are summa- 
rized in Conclusion. All technical details are placed in 
Appendices A - E. 



II. DOUBLE QCB 

A. Basic notions 

Double QCB is a 2D periodic grid, which is formed 
by two periodically crossed arrays of AD quantum wires. 
In experimentally realizable setupsH these are cross- 
structures of suspended single-wall carbon nanotubes 
placed in two parallel planes separated by an inter-plane 
distance d. However, some generic properties of QCB 
may be described in assumption that QCB is a gen- 
uine 2D system. We assume that all wires of j-th array, 
j = 1,2, are identical. They have the same length Lj, 
Fermi velocity vj and Luttinger parameter gj . The arrays 
are oriented along the unit vectors ei^ with an angle tp 
between them. The periods of a crossbars along these 
directions are a\ and 02, and corresponding basic vectors 
are a, = ctjej (Fig.|j). 




FIG. 1. 2D crossbars formed by two interacting arrays of 
parallel quantum wires. Here ei , e2 are the unit vectors of 
the superlattice, a\,a2 are the superlattice periods and d is 
the vertical interarray distance 

The interaction between the excitations in different 
wires is assumed to be concentrated around the cross- 
ing points with coordinates niax + 7T-2&2 = {n\a\, 71202). 



2 



The integers rij enumerate the wires within the j-th ar- 
ray. Such interaction imposes a superperiodicity on the 
energy spectrum of initially one dimensional quantum 
wires, and the eigenstates of this superlattice are charac- 
terized by a 2D quasimomentum q = qigi + (72 g2 = 
(51,52)- Here gi,2 are the unit vectors of the recip- 
rocal superlattice satisfying the standard orthogonal- 
ity relations (e^ ■ gj) — Sij. The corresponding ba- 
sic vectors of the reciprocal superlattice have the form 
™iQigi + miQi&i-, where Qj = 2n/aj and mi. 2 are in- 
tegers. 

However the crossbars kinematics differs from that of 
a standard 2D periodic system. In conventional 2D sys- 
tems, forbidden states in the inverse space arise due to 
Bragg diffraction in a 2D periodic potential, whereas the 
whole plane is allowed for wave propagation in real space, 
at least until the periodic potential is weak enough. A 
Brillouin zone is bounded by the Bragg lines. It coincides 
with a Wigner-Seitz cell of reciprocal lattice. In sharply 
anisotropic QCB most of the real space is forbidden for 
electron and plasmon propagation. The Bragg condi- 
tions for the wave vectors are modulated by a periodic 
potential unlike those in conventional 2D plane. These 
conditions are essentially one-dimensional. Correspond- 
ing BZ is not a Wigner-Seitz cell of a reciprocal lattice 
but the elementary cell containing a site in its center. 

Indeed, the excitation motion in QCB is one- 
dimensional in major part of the 2D plane. The 
anisotropy in real space imposes restrictions on the possi- 
ble values of 2D coordinates x\, x 2 (r = Xie\ +x 2 e2). At 
least one of them, e.g., x 2 (x\) should be an integer multi- 
ple of the corresponding array period a 2 (ai), so that the 
vector r = (x\, n 2 a 2 ) (r = (niai,x 2 )) characterizes the 
point with the ID coordinate x\ (x 2 ) lying at the 712-th 
(ni-th) wire of the first (second) array. As a result, one 
cannot resort to the standard basis of 2D plane waves 
when constructing the eigenstate with a given wave vec- 
tor k. Even in non-interacting arrays of quantum wires 
(empty superlattice) the 2D basis is formed as a super- 
position of two sets of ID waves. The first of them is 
a set of ID excitations propagating along each wire of 
the first array characterized by a unit vector fcigi with a 
phase shift 02^2 between adjacent wires. The second set 
is the similar manifold of excitations propagating along 
the wires of the second array with the wave vector fc 2 g2 
and the phase shift a\k\. The dispersion law of these 
excitations has the form 

w (k) =a; 1 (fc 1 )+w 2 (* 2 ). (2) 

The states of equal energy obtained by means of this 
procedure form straight lines in the 2D reciprocal space. 
For example, the Fermi surface of QCB developed from 
the points ±kp \ >2 for individual quantum wire consists of 
two sets of lines |fci.2| = fc.Fi, 2- Respectively, the Fermi 
sea is not a circle with radius Uf like in the case of free 



2D $as, but a cross in the k plane bounded by these four 
linesE3 (see Fig. ||). Finally, the Bragg conditions read 

w i(^i) — omegai(ki + uiiQi) 

+ ^2(^2) - ^2(^2 + m 2 Q 2 ) = 0. 

and the lines k\ = 0, \k 2 \ = Q 2 /2, and \ki\ — Qi/2, 
k 2 = 0, satisfying these conditions, form a 2D BZ of 
double QCB. 



k 2 




FIG. 2. Fermi surface of 2D metallic quantum bars in the 
absence of charge transfer between wires, gi , g2 are the unit 
vectors of the reciprocal superlattice 

Due to the inter-wire interaction, the excitations 
of QCB (see Figs|^ below) acquire genuine two- 
dimensionality characterized by the quasimomentum q = 
(qi,q 2 ). However, in case of weak interaction the 2D 
waves constructed from the ID plane waves in accordance 
with the above procedure form an appropriate basis for 
the description of elementary excitations in QCB in close 
analogy with the nearly free electron approximation in 
conventional crystalline lattices. It is easily foreknown 
that a weak inter-wire interaction does not completely 
destroy the above quasimomentum classification of eigen- 
states, and the 2D reconstruction of the spectrum may 
be described in terms of wave mixing similarly to the 
standard Bragg diffraction in a weak periodic potential. 
Moreover, the classification of eigenstates of empty su- 
perlattice may be effectively used for the classification of 
energy bands in a real QCB superlattice where the su- 
perperiodicity is imposed by interaction. 

Complete kinematics of an empty superchain (wave 
functions, dispersion laws, relations between quasipar- 
ticle second quantization operators) is developed in Ap- 
pendix A. In terms of these ID Bloch functions (see Eqs. 
(p|), (j44|) of Appendix A) we construct the 2D basis of 
Bloch functions for an empty superlattice 

* SjS ',q(r) = "01,8,81 (xi)ip2,s>, q2 (x 2 ). (3) 
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Here s, s' = 1, 2, . . . , are the band numbers, and the 2D 
quasimomentum q = (91,92) belongs to the first BZ, 
\qj\ < Qj/2. The corresponding eigenfrequencies are 



Here 



uj SS '{q) = wi,s(q) +W2, s '(q)- 



Wj>(q) = Wj>(9j), 



(4) 



and ujj tS (qj) is defined by Eqs. ( pJ5[ ) below. We will use 
this basis in the next subsection when constructing the 
excitation spectrum of QCB within the reduced band 
scheme. 



V(ri — r 2 ) is a short-range interwire interaction. Phys- 
ically, it represents the Coulomb interaction between 
charge fluctuations 

ec( Xj ~ njaj ), C(0 = C(-0, C(0) = 1, (6) 



around the crossing point (nidi, 71202)- The size of the 
fluctuation on the wire of the j'-th array, is determined 
by the screening radius rj within the wire. One may 
neglect the inter-wire tunneling and restrict oneself by 
the capacitive interaction only, provided the vertical dis- 
tance between the wires d is substantially larger than the 
screening radiuses rj. Therefore the interaction has the 
form, 



B. Hamiltonian 

When turning to description of interaction in a QCB, 
one should refer to a real geometry of crossbars, and rec- 
ollect the important fact that the equilibrium distance 
between two arrays is finite and large enough to supress 
direct electron tunneling citeRueckes. We neglect also 
the elastic and van der-Waals components of intertac- 
tion between real nanotubes, because these interactions 
are not involved in formulation of collective excitations 
in QCB. Then the full Hamiltonian of the QCB is 



H — Hi + H 2 + Hi n t , 



(5) 



where Hj describes the II? boson field in the j-th array 

Li/2 



niUi, x 2 ) 



+ — {d Xl 9i (xi,n 2 a 2 )) 2 )■, 
9i 



L2/2 

H 2 = -yY1 j dx 2 js27rf ( 

+ — {dx 2 2 {niai : x 2 )f \, 

92 



and OjiTtj are the conventional canonically conjugated 
boson fields (see, e.g., Ref. ||). 

The interwire interaction results from a short-range con- 
tact capacitive coupling in the crosses of bars, 

Hi n t = E J dxidx 2 V(xi - niai,n 2 a 2 - x 2 ) 

ni ,n 2 

x pi(xi, n 2 a 2 )p 2 (niai, x 2 ), 
where the integration is restricted by the area —Lj/2 < 



w 2 \n ' r 2 
where the function $(^1,^2) is 

Ci(a)c 2 (6) 



(7) 



|£12|_ 

d 2 



1*12 = ri^ex - r 2 £_ 2 e 2 . 

It is seen from these equations that $(£1, £2) vanishes for 
I6-.2I > 1 and is normalized by condition $(0,0) = 1. 
The effective coupling strength is 



V 



2e 2 



(8) 



In terms of boson field operators Qi, the interaction is 
written as 



H int = V E J dx i dx 2® 



xi - niai n 2 a 2 - x 2 



ni ,n 2 

x d Xl &i(xi, 12^)8x^(121^, x 2 ) 



r 2 



In the quasimomentum representation (^|) the full 
Hamiltonian (^) acquires the form, 
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^-E E^wLa 

jj' = l s,s',q 



(9) 



where \/vg/vjgjOj S q and \J Vjgj jvg-Kj S<A are the Fourier 
components of the boson fields 9j and itj , and effective 
velocity and coupling are v = ^/viv 2 , g = y / gig 2 respec- 
tively. 



< Lj/2. Here Pi(r) are density operators, and by: 



The matrix elements for interwire coupling are given 
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Wjsfs'q = Mjs(qj)Uj>a>(qj>) [Sjj'S ga > + fasj's'q (1 ~ %')] ■ 

Here 

wj.fe) = «i([|]0i + (-i)* _1 |«7l). (io) 

are eigenfrequencies of the "unperturbed" ID mode (see 
Eq. ( |45[ ) of Appendix A), pertaining to an array j, band 
s and quasimomentum q = qjgj. The coefficients 

01s2s'q = 0(-f) S+S Signal g 2 )$l S 2s'q, 



hva 



are proportional to the dimensionless Fourier component 
of the interaction strengths 



x < s , gi (n6)w2, s /, g2 (r 2 6) = $2 S 'lsq- 

The Hamiltonian (||) describes a system of coupled 
harmonic oscillators, which can be exactly diagonalized 
with the help of a certain canonical linear transforma- 
tion (note that it is already diagonal with respect to the 
quasimomentum q). The diagonalization procedure is, 
nevertheless, rather cumbersome due to the mixing of 
states belonging to different bands and arrays. However, 
it will be shown below that provided d 3> rip, a separa- 
ble potential approximation is applicable, that shortens 
calculations noticeably. 



C. Approximations 

As it was already mentioned, we consider the rarefied 
QCB with short range capacitive interaction. In the case 
of QCB formed by nanotubes, this is a Coulomb interac- 
tion screened at a distance of the order of the nanotube 
radiusEll Ro , therefore ro ~ Ro ■ The minimal radius of a 
single- wall carbon nanotube is about Ro = 0.35 -r OAnm 
(see Ref. ||^). The intertube vertical distance d in ar- 
tificially produced nanotube networks is estimated as 
d w 2nm (see Ref. ||) Therefore the ratio r 2 /d? w 0.04 is 
really small and the dimensionless interaction 3>(£i,£2) 
(0) in the main approximation is separable 



(11) 



It should be noted that the interaction in this form is an 
even function of its arguments, and the odd correction to 
the $o is of order rjj/d 2 , whereas $o is of order of 1. 

To diagonalize the Hamiltonian (|9|), one should solve 
the system of equations of motion for the field operators. 
Generalized coordinates 9 satisfy the equations 



[wf s (gi) - lo 2 } 6isq + Vs(f>u(qi)oJu(qi) 

X — ^ <p2s' {qi)u2 s > (<72)#2s'q = 0, 



s = l,2,. 



(12) 



and the similar equations obtained by permutation 1 «-> 
2. Here 

cj> js {q) = (-l) s sign(q) / ^(O^S^oO, (13) 



Bloch amplitudes Uj sq (ro£,) are defined by Eqs. ( |44| ) of 
Appendix A, and 



ro) \ hv J 



(14) 



Due to separability of the interaction, equations of mo- 
tion (|l^) can be solved exactly. Corresponding square 
eigenfrequencies are determined by the characteristic 
equation 

F lqi (u 2 )F 2q2 (io 2 ) = - £ , (15) 

where 

Ojtf u%(q)-u> 2 

The function Fj q (cu 2 ) has a set of poles at lo 2 = ui 2 s (q), 
s = 1,2,3,... . For squared frequency smaller than all 
squared initial eigenfrequencies ui 2 s (q), i.e. within the 
interval [0,0^], this is a positive and growing function. 
Its minimal value Fj on the interval is reached at lo 2 = 0, 
and it does not depend on quasimomentum q 

Fj q (o) = = I #c*(0 = p i (1 7 ) 

a J s J 

(here Eqs. ( |l6| ) and ( fl3| ) are used). If parameter e is 
smaller than its critical value 



1 



F1F2 ' 



(18) 



then all solutions lo 2 of the characteristic equation are 
positive. When e increases, the lowest QCB mode soft- 
ens and its square frequency vanishes in a whole BZ 
at e — e c . For exponential charge density distribution 
C(0 = exp(— 1£|), one obtaines e c w 1. 

In our model the dimensionless interaction e in Eq. ( |l4| ) 
can be written as 



(ZRoge 2 

\ d hv 



(19) 



For nanotube QCB, the first factor within parentheses is 
about 0.35. The second one which is nothing but the 
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corresponding QCB "fine structure" constant, can be 
estimated as 0.9 (we used the values of g = 1/3 and 
v = 8 x 10 7 cm/sec, see Ref. |3|). Therefore e approxi- 
mately equals 0.1, so this parameter is really small. Thus 
the considere d sy stem is stable, its spectrum is described 
by Eqs.(15), (|16| ) with a small parameter e. 



corresponding BZ and is also a square (see Fig.g). Reso- 
nant lines are the diagonals of BZ. 



The general Eq. (|l5| ) reduces in infrared limit q, oj — > 
to an equation describing the spectrum of two coupled 
sliding phases, i.e. 1 : 1 arrays in accordance with classi- 
fication , offered in Ref. ^. Equation (3.13) of this paper 
is a long wave limit of our equation ( J47| ) derived in Ap- 
pendix B. Therefore the general analysis of stability of 
the LL fixed point is appicable in our approach. 




D. Spectrum 



Due to the smallness of interaction, the systematics 
of unperturbed levels and states is grossly conserved, at 
least in the low energy region corresponding to the first 
few energy bands. This means that perturbed eigenstates 
could be described by the same quantum numbers (ar- 
ray number, band number and quasimomentum) as the 
unperturbed ones. Such a description fails in two spe- 
cific regions of reciprocal space. The first of them is the 
vicinity of lines qj = nQj/2 with n integer. Indeed, as it 
follows from the equations of motion (^2|), around these 
lines the interband mixing is significant. These lines with 
n = ±1 include BZ boundaries. Because of this BZ which 
is, generally speaking, non relevant, and in this subsec- 
tion we refer mostly to BZ. 



-Q/2 

FIG. 3. Two dimensional BZ of square QCB. 

In the major part of the BZ, for quasimomenta q ly- 
ing far from the diagonals, each eigenstate mostly con- 
serves its initial systematics, i.e. belongs to a given array, 
and mostly depends on a given quasimomentum compo- 
nent. Corresponding dispersion laws remain linear being 
slightly modified near the BZ boundaries only. The main 
change is therefore the renormalization of the plasmon 
velocity. 



The second region is the vicinity of the lines where the 
resonance conditions are fulfilled 



J ls 



(20) 



Here inter-array mixing within the same energy band 
(s = s') or between neighboring bands (s ^ s') is sig- 
nificant. In what follows we will pay attention first of 
all to these two regions because in the rest of the BZ the 
initial systematics of the energy spectrum can be success- 
fully used. 

Equations ([l2]) , (|l5|) , describing the wave functions and 
the dispersion laws are analysed in Appendix B. We de- 
scribe below some of these dispersion qurves for two types 
of QCB basing on this analysis. 
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FIG. 4. The energy spectrum of QCB (solid lines) and non- 
interacting arrays (dashed lines) for quasimomenta at the lines 
OA, FC, and OC of BZ 



1. Square QCB 

We start with the simplest case of square QCB formed 
by identical wires. This means that all parameters 
(wire length, space period, Fermi velocity, LL parame- 
ter, screening radius) are the same for both arrays. The 



In the left part of Fig.[| we displayed dispersion curves 
corresponding to quasimomenta belonging to a generic 
OA line in BZ. In what follows we use (J, s) notations for 
the unperturbed boson propagating along the j-th array 
in the s-th band. Then the lowest curve in this part of 
Fig.|| is, in fact, the slightly renormalized dispersion of 
a (2, 1) boson. The middle curve describes (1, 1) boson, 
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and the upper curve is the dispersion of a (1,2) boson. 
The fourth frequency corresponding to a (2, 2) boson, is 
far above and is not displayed in the figure. It is seen 
that the dispersion remains linear along the whole line 
OA except a nearest vicinity of the BZ boundary (point 
A in Fig§). 

Dispersion curves corresponding to quasi momenta ly- 
ing at the BZ boundary q\ = Q/2, < q 2 < Q/2 (line 
FC in Fig.||) are displayed in the central part in Fig.[|). 
The characteristic feature of this boundary is the intra- 
band degeneracy in one of two arrays. Indeed, in zero 
approximation, two modes (l,s), s — 1,2, propagating 
along the first array are degenerate with unperturbed 
frequency u) = 0.5. The interaction lifts the degeneracy. 
This interaction occurs to be repulsive at the BZ bound- 
aries. As a result the lowest of two middle curves in Fig.|] 
corresponds to (1, u) boson, and upper of them describes 
(1, g) boson. Here the indices g,u denote a boson parity 
with respect to the transposition of the band numbers. 
Note that (l,g) boson exactly conserves its unperturbed 
frequency ui = 0.5. The latter fact is related to the square 
symmetry of the QCB. 

Two others curves correspond to almost non pertubed 
bosons of the second array. The lowest curve describes 
the dispersion of the (2,1) wave. Its counterpart in the 
second band (2, 2) is described by the highest curve in 
the figure. Their dispersion laws are nearly linear, and 
deviations from linearity are observed only near the cor- 
ner of the BZ (point C in Fig.g). 

Consider now dispersion relations of modes with quasi- 
momenta on the diagonal OC of BZ and start with q not 
too close to the BZ corner C {q\ = q% = Q/2). This di- 
agonal is actually one of the resonance lines. Two modes 
in the first band coressponding to different arrays are 
strongly mixed. They mostly have a definite j-parity 
with respect to transposition of array numbers j = 1,2. 
Interaction between these modes occurs to be attractive 
(repulsive) for qiq 2 > (qiq 2 < 0). Therefore the odd 
modes (u, s), at the BZ diagonal OC s = 1, 2, correspond 
to lower frequencies and the even modes (g, s) correspond 
to higher ones. The corresponding dispersion curves are 
displayed in the right part of FigJJ. 

At the BZ corner qi = q 2 = Q/2 (point C in Fig.|) all 
four initial modes j, s = 1 , 2 are degenerate in the lowest 
approximation. This four-fold degeneracy results from 
the square symmetry of BZ (the resonant lines are diag- 
onals of the Z). Weak inter- wire interaction partially lifts 
the degeneracy, however the split modes have a definite 
s-parity with respect to transposition of band numbers 
s = 1,2. The lowest frequency corresponds mostly to 
(g, u) boson, symmetric with respect to transposition of 
array numbers, but antisymmetric with respect to the 
transposition of band numbers. The upper curve de- 
scribes a (u, u) boson with odd both j-parity and s-parity. 



The two middle modes with even band parity, (g,g) and 
(u,g) bosons, remain degenerate and their frequencies 
conserve the unperturbed value w = 0.5. This also re- 
sults from the square symmetry of QCB (m). 
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FIG. 5. Lines of equal frequency of the lowest mode for 
QCB (solid lines) and for noninteracting arrays (dashed lines). 
The lines 1,2,3 correspond to the frequencies U\ = 0.1, 
lo 2 = 0.25, uj s = 0.4 

All these results show that the quantum states of 
the 2D QCB conserve the quasi ID character of the 
Luttinger-like liquid in major part of momentum space, 
and that 2D effects can be successfully calculated within 
the framework of perturbation theory. However, bosons 
with quasimomenta close to the resonant line (diagonal 
OC) of the BZ are strongly mixed bare ID bosons. These 
excitations are essentially two-dimensional, and therefore 
the lines of equal energy in this part of the BZ are mod- 
ified by the 2D interaction (see Fig||). It is clearly seen 
that deviations from linearity occur only in a small part 
of the BZ. The crossover from LL to FL behavior around 
isolated points of the BZ due to a single-particle hy- 
bridization (tunneling) for Fermi excitations was noticed 
in Refs. [l7],[l8|, where a mesh of horizontal and vertical 
stripes in superconducting cuprates was studied. 



2. Tilted QCB 

Now we consider the spectrum of a generic double 
QCB. In this case all parameters (wire length, space pe- 
riod, Fermi velocity, LL parameter, screening radius) de- 
pend, generally speaking, on the array index j. In what 
follows we refer to such a QCB as a tilted QCB. Now the 
resonance condition (po| ) is fulfilled not at the BZ diag- 
onal but at the resonant polygonal line. Its part ODE, 
lying in the first quarter of the BZ, is displayed in Fig.|6| 
(all figures of this subsection correspond to specific val- 
ues V2Q2 = 1, V1Q1 = 1.4). This results in qualitative 
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changes of the spectrum that are related first of all to 
the appearance of two points D and E of the three-fold 
degeneracy for a titled QCB (FigJ|) instead of a single 
point C of four- fold degeneracy for a square QCB (FigJ^). 



q 2 




FIG. 6. BZ of a titled QCB 

We start with the resonant line ODE (Fig.0). The dis- 
persion curves at its OD part and the symmetry proper- 
ties of the corresponding eigenstates are similar to those 
at the OC resonant line for the square QCB (Fig|j|). The 
only difference is that instead of the four-fold degeneracy 
at the BZ corner C of the square QCB, there is a three- 
fold degeneracy at the point D lying at the BZ boundary. 
A completely new situation takes place at the DE line, 
where two other modes (1,1) and (2,2), corresponding 
to different arrays and different bands, are degenerate. 
The interaction lifts this degeneracy and the two middle 
lines in Fig.0 describe even (<?) and odd (u) combina- 
tions of these modes. The even mode corresponds to the 
lowest frequency and the odd mode corresponds to the 
higher one. At the point E one meets another type of 
a three-fold degeneracy described in more detail in the 
next paragraph.. 

Dispersion curves corresponding to quasi momenta ly- 
ing at the BZ boundary q x = Qi/2, < q 2 < Q 2 /2 (FC 
line in Fig.|) and q 2 = Q 2 /2, < q 1 < Qi/2 (CF' line in 
Fig.^), are displayed in Fig.^. The lowest and the high- 
est curves in the FE part of the latter figure, describe 
two waves propagating along the second array. They are 
nearly linear, and deviations from linearity are observed 
only near the point E where the interaction has a reso- 
nant character. Two modes propagating along the first 
array, in zero approximation, are degenerate with an un- 
perturbed frequency u> = 0.7. The interaction lifts the 
degeneracy. The lowest of the two middle curves corre- 
sponds to (l,u) boson, and the upper of one describes 



(1, g) boson. Note that (l,g) boson conserves its unper- 
turbed frequency uj = 0.7. The latter fact is related to 
the symmetry Ci(£) = °f the separable interac- 
tion (ph. At the point E, the two modes propagating 
along the first array and the mode propagating along the 
second array in the second band are degenerate. Interac- 
tions lifts the degeneracy, and, as a result, the (1, u) and 
(2, 2) waves are strongly mixed and the eigenmodes are 
their even (highest frequency) and odd (lowest frequency) 
combinations, and the (l,g) mode (middle level). 
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FIG. 7. The energy spectrum of a tilted QCB (solid lines) 
and noninteracting arrays (dashed lines) for quasimomenta on 
the resonant line of the BZ (line ODE in Fig.^J) 
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FIG. 8. Energy spectrum of a tilted QCB (solid lines) and 
noninteracting arrays (dashed lines) for quasimomenta at the 
BZ boundary (line FCF' in the FigJ|) 

There are two separate degeneracies within each array 
at the corner C of a titled QCB BZ. Both of them are 
related to interband mixing conserving array index. The 
spectral behavior along the CF 1 boundary of the BZ is 
similar to that considered above but in the vicinity of the 
point D of three-fold degeneracy. Here, two modes prop- 
agating along the second array in the separable potential 
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approximation (|llj) remain degenerate. This degeneracy 
is lifted only if deviation from separability is accounted 
for. 
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FIG. 9. The energy spectrum of a titled QCB (solid lines) 
and noninteracting arrays (dashed lines) for quasimomenta on 
the BZ diagonal (line OC in Fig.|) 




FIG. 10. Lines of equal frequency for a tilted QCB (solid 
lines) and noninteracting arrays (dashed lines). Lines 1,2,3 
correspond to frequencies u>i — 0.1, u>2 ~ 0.25, UJ3 = 0.45. 



The diagonal OC of a tilted QCB BZ represents a new 
type of generic line, that crosses a resonant line (FigJ|). 
Here the spectrum mostly conserves its initial systemat- 
ics, i.e. belongs to a given array, and mostly depends 
on a given quasimomentum component. However, at the 
crossing point B, the modes (1,1) and (2,2), correspond- 
ing to both different arrays and bands, become degener- 
ate (two middle dashed lines in FigJ|). Interaction be- 
tween the wires lifts the degeneracy. The eigenstates of 



QCB have a definite parity with respect to transposition 
of these two modes. The lowest and upper of two middle 
lines corresponds to even (g) and odd (u) mode, respec- 
tively. 

Like in square QCB, bosons with quasimomenta close 
to the resonant lines are strongly mixed bare ID bosons. 
These excitations are essentially two-dimensional, and 
therefore lines of equal energy in the vicinity of the 
resonant lines are modified by the 2D interaction (see 
Figs. [To| , |TT| ) . Deviations from ID behaviour occur only 
in this small part of the BZ. For uj < 0.5«2<52 the lines 
of equal energy within BZ consist of closed line around 
the BZ center and four open lines (within the extendcnd 
bands scheme these lines are certainly closed) around the 
BZ corners (lines 1, 2, 3 in Figjlp|). At the line OD in 
BZ, the modes of QCB are strongly coupled bare bosons 
propagating along both arrays in the first band. 




FIG. 11. Lines of equal frequency for a tilted QCB (solid 
lines) and noninteracting arrays (dashed lines) . Lines 4, 5 
in the lower panel correspond to frequencies U14 = 0.55, 
us = 0.65 



For 0.5u 2 <32 < lo < 0.5i>i<3i (lines 4, 5 in Fig.^TJ) the 
topology of lines of equal energy is modified. In this case 
lines of equal energy within the BZ consist of four open 
lines. The splitting of lines at the direction DE corre- 
sponds to strong coupling of modes propagating along 
the first array in the first band with those propagating 
along the second array in the second band. 



E. Correlations and Observables 

The structure of the energy spectrum analyzed above 
predetermines optical and transport properties of QCB. 
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We consider here three types of correlation functions 
manifesting dimensional crossover in QCB. 



1. Optical Absorption 

We start with oc conductivity 

a r y (q, uj) = o-' jr (q, uj) + ia'^, (q, to). (21) 

The real part a'^, (q, uj) determines an optical absorption. 
The spectral properties of oc conductivity are given by a 
current-current correlator 

oc 

vjj, (q, uj) = i f dte™ ( [j jlq (t), Jj, lq (0)] ) . (22) 



Here Jj sq = \2vg-Kj stl is a current operator for the j'-th 
array (we restrict ourselves to the first band, for the sake 
of simplicity) . 

The current-current correlator for non-interacting 
wires is reduced to the conventional LL expression^, 



with metallic-like peak 



,(q, uj > 0) = nvg8(uj - u) jltl )Sj 



(23) 



For QCB this correlator is calculated in Appendix C. Its 
analysis leads to the following results. 

The longitudinal absorption 

<7£i(q,w) oc (1 - 4> 2 1ci )5(uj - u>i q ) + 4> 2 ltl 8{u) - <Z> 2q ) (24) 

contains well pronounced peak on the modified first array 
frequency and weak peak at the second array frequency 
(the parameter <^>i q , defined by Eq. (|58|)of Appendix B, 
is small). The modified frequencies a>i q and uj 2<i coincide 
with the eigenfrequencies w+i q and a->_2 q respectively, if 
cji q > uj2q- In the opposite case the signs +, — should be 
changede to the opposite ones. 

The transverse absorption component contains two 
weak peaks 



a' 12 (q,uj) oc 0i q [S(uj - u>i q ) + S(uj - w 2q )] 



(25) 



At the resonant line, the results change drastically. 
Both longitudinal and transverse components of the op- 
tical absorption contain two well pronounced peaks cor- 
responding to slightly split modified frequencies 

cr' n (q,uj) oc i [S(uj - Wi q ) + S(uj - w 2q )] ■ (26) 



2. Space Perturbation 

One of the main effects specific for QCB is the ap- 
pearance of non-zero transverse momentum-momentum 
correlation function. In space-time coordinates (x, t) its 
representation reads, 



G12 (x, t) = i ([tti (x! , 0; t), 7T 2 (0, x 2 ; 0)]) 



(27) 




FIG. 12. The transverse correlation function Gi2(xi, X2',t) 
for ro = 1 and «t = 10 

This function describes the momentum response at the 
point (0,a;2) of the second array for time t caused by 
an initial (t = 0) perturbation localized in coordinate 
space at the point {xi,0) of the first array. Standard 
calculations similar to those described above, lead to the 
following expression, 



Gi 2 (x;t) 



dk 1 dk24>i(ki)(f>2(k 2 )kxk2 



x sin(fcixi) sin^a^) 



^fe sin(i>2fc2£) — v iki sin(wifcit) 
v\k\ - vfkf 



where (j)j(k) is the form-factor ( p^[ ) written in the ex- 
tended BZ. This correlator is shown in Fig.[l2]. It is 
mostly localized at the line determined by the obvious 
kinematic condition 



Fil 



F2| 

v 2 



= t. 



The time t in the r.h.s. is thea total time of plasmon 
propagation from the starting point {x\ , 0) to the final 
point (0, x 2 ) or vice versa, along any of the shortest ways 
compatible with a restricted geometry of the 2D grid. 
The finiteness of the interaction radius slightly spreads 
this peak and modifies its profile. 
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3. Rabi Oscillations 



A. Notions and Hamiltonian 



Further manifestation of the 2D character of QCB sys- 
tem is related to the possibility of periodic energy transfer 
between the two arrays. Consider an initial perturbation 
which excites a plane wave with amplitude 9q within the 
first array in the system of now-interacting arrays, 

9i{xi,n 2 a 2 ;t) = o sm(9ixi -I- q 2 n 2 a 2 - i'i\qi\t). 

If the wave vector q, satisfying the condition |q| << 
Qi^/2, is not close to the resonant line of the BZ, weak 
interarray interaction <p = er /a slightly changes the 

91 component and leads to the appearance of a small 

9 2 ~ 4> component. But for q lying on the resonant line 
(«i = = u q ), both components within the main 
approximation have the same order of magnitude, 



7 ( 1 . lh><i 2 : I ) = 9 n cos ( ^iqWqt 



x sin(gia;i + q 2 n 2 a 2 - w q i), 



y 2 {niai,x 2 ;t) = 6»o sin ( -0i q w q i 



x cos(qiniai + q 2 x 2 - uj^t). 




FIG. 13. Periodic energy exchange between arrays (Rabi 
oscillations) 

This corresponds to 2D propagation of a plane wave 
with wave vector q, modulated by a "slow" frequency 
~ (fxjj. As a result, beating arises due to periodic energy 
transfer from one array to another during a long period 
T ~ ((jyuj)^ 1 (see Fig.|l3|) . These peculiar "Rabi oscilla- 
tions" may be considered as one of the fingerprints of the 
physics exposed in QCB systems. 



Triple quantum bars is 2D periodic grid with m = 3, 
formed by three periodically crossed arrays j = 1, 2, 3 of 
II? quantum wires. In fact these arrays are placed on 
three planes parallel to XY plane and separated by an 
inter-plane distances d. The upper and the lower arrays 
correspond to j = 1,2, while the middle array has num- 
ber j — 3. All wires in all arrays are identical. They 
have the same length L, Fermi velocity v and Luttingcr 
parameter g. The arrays are oriented along the 2D unit 
vectors 



ei 




e 2 



(1,0), 



e 3 = e 2 - ei. 



(28) 



The periods of QCB along these directions are equal, 
Oj = a, so we deal with a regular triangular lattice. In 
what follows we choose ai 2 — aei, 2 as the basic vectors 
of a superlattice (see FigQ. 

The wires within the j'-th array are enumerated with 
the integers nj. Define 2D coordinates along the nj-th 
wire Tj as Yj — Xjej + njae^ for upper and lower arrays 
(j = 1,2) and r3 = 2:363 + n^aei for the middle array. 
Here Xj are ID continuous coorinates along the wire. 
The system of three non-interacting arrays is described 
by the Hamiltonian 



Ho — Hi + H 2 + Hs 



(29) 



where 



Hi = -y ^ J dx i [9^1 ( x i e i + niae 3 ) 



1 2 

+ - (d Xl 9 1 (x 1 e 1 + mae 3 )) 
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Ho = 



H. = 



^ [ dx 2 [g-Kl(x 2 e 2 + n 2 ae 3 ) 

1 2 " 

+ - (d X2 9 2 {x 2 e 2 + n 2 ae 3 )) , 
9 

J dx3 [ ff7I 3 ( :E3e3 + "saei) 



1 2 

+ - (d X3 9 3 (x 3 e 3 + n 3 aei)) 



(30) 



(31) 



(32) 



and ttj and d Xj 9j are canonically conjugated fields de- 
scribing LL within the j-th array. 



III. TRIPLE QCB 
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Interaction between the excitations in different wires 
of adjacent arrays is concentrated near the cross- 
ing points with coordinates UjBj + Uj/a,/. It is actually 
Coulomb interaction screened on a distance ro along each 
wire which is described by Hamiltonian 



Hint — Hi 3 + H 23 , 



(33) 



where 



if = £ / ***** 



§2 

FIG. 15. Elementary cell BIJL of the reciprocal lattice 
and the BZ hexagon of the triple QCB 

To study the energy spectrum and the eigenstates of 
the total Hamiltonian 

H = H + H int , (37) 

we define the Fourier components of the held operators 



xi — n 3 a x 3 — n\a 

ei e 3 



ro 



ro 



x d Xl Q\{x\&\ + niae 3 )d X3 9 3 (n 3 ae 1 + x 3 e 3 ) 



(34) 



if = £ / ***** 



X2 — n 3 a x 3 — n 2 a 

e 2 e 3 



ro 



ro 



x d X2 02(x 2 e 2 + n 2 ae 3 )d X3 9 3 (n 3 ae 2 + x 3 e 3 ). (35) 

Here the effective coupling strength Vq is defined by 
Eq.(|J), the dimensionless interaction <E> is separable 

*& ej + £ 3 e 3 ) = C(0)C(6), 3 = 1, 2, (36) 

and is dimensionless charge fluctuation in the j-th 
wire (see Eq. (ph). 

Such interaction imposes a superperiodicity on the 
energy spectrum of initially one dimensional quantum 
wires, and the eigenstates of this superlattice are char- 
acterized by a 2D quasimomcntum q = (?igi + <Z2g2 = 
(qi,q 2 ). Here gi^ are the unit vectors of the reciprocal 
superlattice satisfying the standard orthogonality rela- 



0i(xiei + niae 3 ) = 
{NL)- 1 / 2 Y.Qis^ qixl+q3n ^u s , qi {x 1 ) 1 

Q 2 (x 2 e 2 + n 2 ae 3 ) = 

(7Vi)- 1 /2^0 2sqe »(92^+ 9 3n 2 a) Us(j2(x2); 

d3{x 3 e 3 + n 3 aei) = 
(NL)- 1 ' 2 Y, 0^e l ^+^ n ^u sm {x 3 ). 



(38) 
(39) 
(40) 



Here 



q = qiei + q 2 e 2 , q 3 = q 2 - qi, 



and N = L/a is the dimensionless length of a wire. In the 
q representation, the Hamiltonians H.j (Eqs. (|30|)-(32)) 
and Hj 3 (Eqs. (pi]) , (^)) can be written as 



hvg 



J = 1,2,3, 



tions (ej • gj) 



J 



1,2. The corresponding ba- 



H 



sic vectors of the reciprocal superlattice have the form 
Q(migi + m.2g2, where Q = 2n/a and mi^ are integers. 
In Fig.p^| elementary cell BIJL of the reciprocal lattice is 
displayed together with the hexagon of the Wigner-Seitz 
cell that we choose as the BZ of the triple QCB. 



23 



V r 2 



2v 9 f~ 



£ <f>s{q3)^s'(qj)^ s (q3)^s'(qj) 



8,8' ,q 



where 



K^ 3 s^ + h.c], j = 1,2, 



(q)=v{[^\Q+(-ir 1 \q\). Q = ^ 
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Thus the total Hamiltonian (^) describes a system of 
coupled harmonic oscillators, and can be diagonalized ex- 
actly like in the case of double QCB. 



B. Spectrum 

Separability of the interaction (j3(^) allows one to de- 
rive analytical equations for the spectrum of the total 
Hamiltonian ( p7| ) (see Appendix D). Here we describe 
the behavior of the spectrum and the states along some 
specific lines of the reciprocal space. 

The high symmetry of the triple QCB leads to a num- 
ber of lines where interarray or interband resonant inter- 
action occurs: all lines in Fig. [l5] posess some resonant 
properties. These lines may be classified as follows: 

On the Bragg lines where one of three array wavenum- 
bers qj is a multiple integer of Q/2, there is a str ong 
intraband mixing of modes of the j-th array. In Figjlq, 
these lines are the boundaries of the elementary cell of 
the reciprocal lattice IJLB, axes qi and q2, lines OB 
and EH. In particular, along the lines OA (q2 = 0) and 
OB (q 3 = 0) two modes corresponding to 2-d and 3-d 
bands and to the second (OA) or third (OB) array are 
mixed. Along the line AB (qi — Q/2) the same mixing 
happens between (1,1) and (1,2) modes. Moreover, the 
resonant mixing of different arrays within the same band 
occurs along the medians OA, OB, etc.. There are two 
types of such a resonance. The first of one (e.g. OA line) 
is the resonance between neighboring arrays (q± = —qi,) 
and therefore it is of the main order with respect to inter- 
action. The second one (e.g. OB line) is the resonance 
between remote arrays (gi = qi) and is one order smaller. 



types of two-fold degeneracy simultaneously. For exam- 
ple, at the point M two separate pairs of modes corre- 
sponding to neighboring arrays (2,1), (3,1), and (2,2), 
(3, 2), are degenerate, as well as two modes corresponding 
to the first array, (1,1), (1,2). Finally the BZ hexagon 
vertices form the most interesting group of points where 
the three-fold degeneracy between modes corresponding 
to all three arrays takes place. The typical example of 
such a point is the vertex C where the resonance condi- 
tion qi = —q 2 + Q = qs = Q/3 is satisfied. 




O A M B O 



FIG. 16. Dispersion curves at the OAMBO polygon of BZ 

Almost all these peculiarities of the triple QCB spec- 
trum can be illustrated in Fig.|l6| where the dispersion 
curves along the closed line OABO are displayed. We 
emphasize once more that in the infrared limit lo, q — * 
triple QCB like double QCB preserves the characteristic 
LL properties of the initial arrays. 



The second family consists of resonant lines formed by 
the BZ hexagon boundaries and diagonals. Thus, the di- 
agonal OC realizes a first order resonance between the 
first and the third arrays q\ — qj,, and the BZ boundaries 
HD and AN correspond to the same resonance up to an 
umklapp process (q± = q^ — Q and q\ = q% + Q respec- 
tively). Along the diagonal OD and the BZ boundary 
NC a second order resonance takes place with resonance 
conditions qi — —q± and 52 = — Qi + Q respectively. 

In the reciprocal space of the triple QCB there are 
four different types of crossing points. Two of them in- 
clude the bases of BZ medians (e.g. points A, B, E and 
so on). Here one deals with the four- fold degeneracy of 
the modes corresponding to the first order resonance be- 
tween the neighboring arrays (e.g. point A, u>i jS = 073, s /, 
s,s' = 1,2), or to the second order resonance between 
remote arrays (like point B, u>±^ s = u>2, s >, s,s' = 1,2). 
One more family consists of crossing points of the BZ di- 
agonals and the lines connecting the bases of its medians 
(points M, F, G and so on). Here one deals with three 



C. Observables 

The structure of the energy spectrum analyzed above 
strongly influences optical and transport properties of the 
tripl e QCB. As in the case of double QCB (subsection 
|IIE|) , one expects to observe four peaks of the optical 
absorption near the points A, B, E, H of the four-fold de- 
generacy. Then, specific features of space correlators like 
those considered in HE 2 can be observed. But the most 



pronounced manifestation of a triangular symmetry of 
the triple QCB are its Rabi oscillations. 

Consider the vicinity of the point C of three-fold de- 
generacy mixing all three arrays. Appropriate initial con- 
ditions lead (see Appendix E for details) to the following 
time dependence of the field operators in the coordinate 
origin in real space 



0i(O,O;t) = 9 sin(w t) cos 2 



2e(j)< 



u> t 
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2 (0,0;*) = 6> cos(w *) sin 2 



9 3 (0, 0; *) = 8 sin(w *) cos 



2e<p< 



UJ t 



2s<f> 2 



(41) 



The field operators of all three arrays demonstrate fast 
oscillations with the resonant frequency a->o modulated 
by a slow frequency. It is the same for the two remote 
arrays, and doubled for the intermediate array. These 
beatings are synchronized in a sense that zero intensity 
on the intermediate array always coincides with the same 
intensity on one of the remote arrays. At these moments 
all the energy is concentrated solely within one of the 
remote arrays. These peculiar Rabi oscillations are dis- 
played in Fig.|l7|. 
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FIG. 17. Periodic energy transfer between three arrays at 
the triple resonant point C of the BZ 



and confirmed in the present study, the energy spectrum 
of QCB preserves the characteristic properties of LL at 
\q\,u> — ► 0,. At finite q, u> the density and momentum 
waves in double and triple QCB may have either ID or 
2D character depending on the direction of the wave vec- 
tor. Due to interarray interaction, unperturbed states, 
propagating along arrays are always mixed, and trans- 
verse components of correlation functions do not vanish. 
For quasi-momentum lying on the resonant lines of the 
BZ, such mixing is strong and transverse correlators have 
the same order of magnitude as the longitudinal ones. 
Periodic energy transfer betweem arrays ("Rabi oscilla- 
tions" ) is predicted. 

The crossover from ID to 2D regime may be exper- 
imentally observed. One of the experimental manifes- 
tations, i.e. the crossover from isotropic to anisotropic 
(spatially nonuniform) conductivity was pointed out in 
Ref. 0. The current may be inserted in QCB at a point 
on an array j and extracted from another array i at a 
distance r. Then a temperature dependent length scale 
l(T) arises, so that for r ^ / the resistance is dominated 
by small q and therefore, the current is isotropic. In the 
opposite limit r < I the dependence of the current on 
the points of injection/extraction may be detected. At 
T = the length I becomes infinite, and current can only 
be carried along the wires. These effects are in fact man- 
ifestations of the LL behavior of the QCB in the infrared 
limit. 

To observe the crossover at finite {o>,q}, one should 
find a way of exciting the corresponding plasmon modes. 
Then, scanning the Lo{qi,q2) surfaces, one may in prin- 
ciple detect the crossover from quasi ID to 2D behavior 
in accordance with the properties of the energy spectra 
presented in Sections II and III. Plasmons in QCB may 
be excited either by means of microwave resonators or by 
means of interaction with surface plasmons. In the lat- 
ter case one should prepare the grid on a corresponding 
semiconductor substrate and measure, e.g., the plasmon 
loss spectra. The theory of these plasmon losses will be 
presented in a forthcoming publication. 
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APPENDIX A. EMPTY SUPERCHAIN 

Here we construct eigenfunctions, spectrum, and 
quasi-particle operators for an "empty superchain" - 
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quantum wire in an infinitely weak periodic potential 
with period a. Excitations in an initial wire are de- 
scribed as plane waves L -1 / 2 exp(ikx) with wave num- 
ber k = 2nn/L, with integer n, and dispersion law 
u(k) = v\k\ (the array number is temporarily omitted). 
The following orthogonality relations are valid 



,L/2 
/ 

J-L/2 



ipt( x )^k'{x)dx = 5 k ,k', 
^2^* k {x)ip k {x') = 5 L (x-x'), 

k 

where Sl stands for periodic delta- function 

Sl(x — x') = 6(x — x' — nL). 

n 

"Empty superchain" is characterized by a space pe- 
riod a and corresponding reciprocal lattice wave number 
Q = 2n/a. Each excitation in such a superchain is de- 
scribed by its quasi-wavenumber q and a band number s 
(s = 1,2,...) that are related to the corresponding wave 
number k by the following relation, 



k = q + iQx(-l) s 



sign q. 



(42) 



The corresponding wave function ip Stq (x) has the Bloch- 
type structure, 



V> s>9 (z) = -^=e iqx u^ q (x), 
and satisfies the orthogonality relations 

r L/2 



(43) 



L/2 



^*s,q{x)^8',q'{x)dx = 5 s ^5 Q . q ^, 
^^*s,q{x)lp s , q {x') =S L (X-X'), 



s,q 



where 



$Q;q,q' - ^+™ f « 



,q'- 



Within the first BZ, -Q/2<q< Q/2, Bloch amplitude 
and dispersion law to s have the following form 



u s , q ( x ) = exp|iQa;(-l) s 1 | sign qj , 



(9) = «Q([f]+(-i)'- 1 



(44) 
(45) 



Here square brackets denote an integral part of a number. 
Taking into account that both Bloch amplitude u s , q (x) 
and dispersion law ui s (q) are periodic functions of q with 
period Q, one obtaines general equations for the Bloch 
amplitude 



sin^„ 



x cos[(2s - 1)£„] exp ( -4i£„ — 

4£n = Q(x - na). 
and dispersion law w s {q) 



{vQ)- 1 uj s {q) = 



2s- 1 



2(-l)« 



^7T 2 (2n + l) 2 



COS ■ 



27r(2n + l)q 

~Q ' 



The relations between quasiparticle operators for a free 
wire, Cfe, for momentum k ^ nQ/2 with n integer, and 
those for an empty superchain, C s>9 , for quasimomentum 
q from the first BZ, -Q/2 < q < Q/2, look as 



Q 



k = q+{-l) v Q 



Ck = Cg.qSignk, 

°<{HH 

C s ,q — {— 1) I/ C/ S , 

-2q] 



s + 1 



Q 



where curely brackets denote a fractional part of a num- 
ber. For obtaining these relations we used the folowing 
expression 

r L/2 

^t(x)ips,q(x)dx = 5 s , s ( q) 5 Qiq . k signk, 



L/2 



s(q) = 1 + 



2kl 



Q 



for the transition amplitude (k\s,q). In case when k = 
nQ/2 with n integer, hybridization of the neighboring 
bands should be taken into account. This modifies the 
above relations by the following way 

c nQ/2 = 0(n) [a n C n ,q n + P n C n +l,q n ] 

+9(-n) [p*_ n C- n , qn - a*_ n C- n+1 , qn ] , 



Q 



C s ,q s — a*C s Q/2 + PsC^ s Q/2, 



1,9. 



KCsQ/2 



a s c- 



sQ/2, 



where a, (3 are hybridization coefficients. Correspond- 
ing relations between wave functions follow immediately 
from these formulas. 

To write down any of these formulas for a specific ar- 
ray, one should add the array index j to the wave function 
tp, Bloch amplitude u, coordinate x, quasimomentum q, 
and to the periods a and Q of the superchain in real and 
reciprocal space. 
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APPENDIX B. DOUBLE QCB SPECTRUM 

Here we obtain analytical expressions for dispersion 
laws and wave functions of QCB. For quasimomenta far 
from the BZ boundaries, the energy spectrum of the first 
band can be calculated explicitly. Assuming that u 2 <C 
uj 2 8 (<7j), s = 2,3,4,..., we omit co 2 in all terms in the 



r.h.s. of Eq. (16) except the first one, s = 1. As a result 



the secular equation (17) reads 



n 



=i y^ji^)-^ 



(46) 



where 



a j 

The solutions of this equation have the form: 2 

2 

/lq 



^i q = w 2 (q)+w 2 (q) 



Here 



± ^(q) -J> 2 2 (q)) 2 +4 £ ^?(g 1 H 2 ( g2 ). (47) 



Vq = <<5l('7l)^2(<Z2), 

is the branch number, cDj(q) is determined as 



«i(q) =uM 



l-e^-^gi))^-^)) 



(48) 



for j=l. Expression for wf(q) can be obtained by per- 
mutation 1 <-> 2. Parentheses on the r.h.s. of Eq. fljsj ) de- 
scribe the contributions to Fj from higher bands. There- 
fore Cj 2 (q) is the j-th array frequency renormalized by the 
interaction with higher bands. In principle, contribution 
of higher bands may turn the interaction to be strong. 
However for specific case of carbon nanotubes, one stays 
far from th e cri tical value e c (see estimates at the end of 
subsection II C). Therefore the interaction with higher 



bands is weak almost in all the BZ except its boundaries. 

The resonance line equation modified by interaction 
with higher bands is 



Wi(q) =££(q 



(49) 



Out of this line the branch number is in fact the array 
number and the renormalized frequencies are frequencies 
of a boson propagating along one of the arrays slightly 
modified by interactions with the complementary array. 
In case when u)\{qx) > 0^2(92)1 one obtains 

<iq«w?( gi )(l-ei^(«i)). ( 5 °) 
In the opposite case one should replace indices 1 2, 



Consider the frequency correction in the latter equa- 
tion in more details. The correction term can be approx- 
imately estimated as uj\{qi)S{qi) with 



S( qi ) = eF 2 rf( qi ) = e^4> 2 n (qi) I «Cf (0- (51) 



Due to the short-range character of the interaction, the 
matrix elements <frn(qi) ~ 1 vary slowly with the quasi- 
momentum q\ < Q\. Therefore, the r.h.s. in Eq.(pl|) can 
be roughly estimated as 



S( gi )~3 = 0.1^<l. 
a a 



(52) 



One should also remember that the energy spectrum 
of nanotube remains one-dimensional only for frequen- 
cies smaller than some u> m . Therefore, an external cutoff 
arises at s = ak m where k m ~ tv m /v. As a results one 
gets an estimate 



S{qi) ~ e—k m R . 
a 



(53) 



Hence, one could hope to gain additional power of the 
small interaction radius. However, for nanotubes, k m is 
of the order of 1/Rq (see Refs. |26| , p3| ) and both estimates 
coincide. For quasimomenta close to the BZ center, the 
coefficient S(qi) can be calculated exactly. For exponen- 
tial form of oc exp(— 1£|), one obtains instead of the 
preliminary estimate ([52]), 

5(0) = 0.14—. 

a 

Thus, the correction term in Eq.(|50|) is really small. 

The eigenstates of the system are described by renor- 
malized field operators. Within the first band they have 
the form 



'llq 



1 - ^/?lq ) (UqOnq - Vq8 2 lq) 



E 1 

s=2 



ilsqMq#2sq + 02sqVq#lsq) , 



721q 



1 ~ 2^ 2t J ) ("q^Hq + u q^21q) 



E 1 

s=2 



ksq^qfesq + ^sqUq^lsq) 



(54) 



(55) 



Here the coefficients it q and w q describe mixing between 
the modes with different array indices, within the first 
band, 



\ 



2,/A q 



(56) 
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\ 



'lq 



2./A q - 



(57) 



and 



01q = V£ ( Pq^ll('ll)^2l(q2)- 

The parameters </>i sq , </>2sq, s = 2,3,. 
(l55|) correspond to inter-band mixing 



(58) 



in Eqs. (54) 



5i sq = ye— (f>n{qi)(P2s{q2)- 



U2s{q2) ' 



(59) 



and the coefficients /3i q , /3 2q , take into account correc- 
tions from the higher bands 



lq 



E 

s=2 



5 2s q«q 



(60) 



Expressions for 02s q and /?2q can be obtained by permu- 
tation 1^2. 



Equations (^7|), ( j54| ) and (55) solve the problem of 
QCB energy spectrum away from the BZ boundaries. 
However, due to smalness of the interaction, the general 
expressions (|54j) and ( |55| ) can be simplified. For quasimo- 
menta far from the resonant coupling line, the expressions 
for the renormalized field operators of the first array look 
like 



?llq 



where 



1 



i>ii q = 0ii (91)^21 {q%/ ■■>,-■ 



1 



lq I Cllq 



s=l 



'lsqC2sq, 



(61) 



wii(gi)w 2 i(g 2 ) 



and 



Aq = E^q- 



(62) 



The corresponding formulas for the second array are ob- 
tained by replacing Is — > 2s. 

Another simplification is made for modes with quasi- 
momenta on the resonance line. Consider for simplicity 
a square QCB (in this case BZ coincides with the ele- 
mentary cell of the reciprocal lattice, and the resonance 
line coincides with the BZ diagonal OC in Fig.|J) and as- 
sume that q is not too close to the BZ corner C. The 
initial frequencies of modes belonging to the same band 
coincide, 



Wlsq — W2sq 



y sq ■ 



(63) 



Therefore renormalization strongly mixes the initial vari- 
ables 



?gsq = —7= I 1 — ^/3sq I ($2sq + #lsq) 



—7= (4>s'sq9ls'q ~ 4*ss'q9 



'qC 2 s'qJ, 



?tisq — —7= I 1 — ^/3sq I ~ ClsqJ 



J]] {4>s>sq8ls'ci + 0ss'q#2s'c 



and the corresponding eigenfrequencies are shifted from 
their bare values 



W sq (1 + 01s2 S q) 



W «sq ~ W fq (1 - 01s2sq) • 

For the first band s = 1 these formulas look like 

< lq «wi 2 (5i)(l±^ q ). (64) 

Note that in the resonance case the splitting of the de- 
generate modes is of the order of \fe that essentially ex- 
ceeds the shift of eigenfrequencies in the non-resonant 
case @. 

The interband mixing becomes significant near the BZ 
boundaries. Not very close to the crossing points of these 
boundaries with the resonant lines, this mixing is ac- 
counted for by a standard way. As a result we find that 
the interband hybridization gap for the bosons propagat- 
ing along the first array can be estimated as 



Auj 



12 



vQe—. 
a 



Similar gaps exist near the boundary of the BZ for each 
pair of odd and next even energy bands, as well as for 
each even and next odd band near the lines qi = or 
#2 = 0. The energy gap between the s-th and (s + l)-th 
bands is estimated as 



Auj 



s,s + l 



?e%(s- 1 ). 
a 



For large enough band number s, the interaction is effec- 
tively suppressed, 4>is2s' — > 0, and the gaps vanish. 

The spectral behavior in the vicinity of the crossing 
points of a resonance line and the BZ boundary needs 
more detailed calculations. Nevertheless it can also be 
analyzed in a similar way. The result s of such an alalysis 
are discussed in subsubscction II D 1 . 
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APPENDIX C. AC CONDUCTIVITY 

For interacting wires, where <pj S {qj) ^ 0, the correlator 
( [22| ) may be easily calculated after diagonalization of the 
Hamiltonian ([)]) by means of the transformations ( p34| ) 
and (Ba). As a result, one has: 



■/Hq(<),4lq(0) 

-2wg (jiqW +i i q sin(cj +i i q £) + VqW-,i q sin(w_ i q t)) 



-2ivgUqV q 



x (o;_,i q sin(a;_ i iqt) - w+,i q sin(w +4q i)) . 



where u q and D q are defined in Eqs.(p6|), (57). Then, for 
the optical absorption a' one obtains 

cr^q.a;) = nvg u 2 S (w - £+,i q ) + v q <5 (w - w_,i q ) (65) 



(7^2 (q,w) = 7rugu q z; q <5(w - £-,i q ) - 5 (w - £+,i q ) 



(66) 



<(q,-) = ^ 



2 2 2 2 

«,lq 
o o "T" o 



Iq 



,iq 



-,iq 



Beside the standard pole at zero frequency, the imaginary 
part has poles at the resonance frequencies w+^q, w_,i q , 
and an additional series of high band satellites (omitted 
here). For quasimomenta far from the resonant lines, 
only the first pole is well pronounced while amplitude of 
the second one as well as amplitudes of all other sattelites 
is small. At the resonant lines, amplitudes of both poles 
mentioned above are equal. The corresponding expres- 
sion for (j' 2 2((i, w) can be obtained by replacement 1 <-> 2. 

The transverse component of the imaginary part of the 
ac conductivity has the form: 



<?12 (q> W ) = W q U q 

id 



■lq 



-,iq 



,lq 



It always contains two poles and vanishes for noninter- 
acting wires. For quasimomenta far from the resonance 
lines the transverse component is small while at these 
lines its amplitude is of the order of unity. 



For quasimomentum q away from the resonant coupling 
line, w q « 1 and ~ q for A q > (u q « 1 and u 2 



for A q < 0). Then the longitudinal optical absorption 
( |65| ) (i.e. the absorption within a given set of wires) has 
its main peak at the frequency w+,iq ~ v\q\\ for A q > 
(or w_ : i q w v\qi\ for A q < 0), corresponding to the first 
band of the pertinent array, and an additional weak peak 
at the frequency W-,iq ~ f|<?2|i corresponding to the first 
band of a complementary array. It contains also a set of 
weak peaks at frequencies W2, sq ~ [s/2]vQ (s = 2, 3, . . .) 
corresponding to the contribution from higher bands of 
the complementary array (in Eq.(|65]) these peaks are 
omitted). At the same time, a second observable becomes 
relevant, namely, the transverse optical absorption (66). 
It is proportional to the (small) interaction strength and 
has two peaks at frequencies w+,iq and u>_.i q in the first 
bands of both sets of wires. 

If the quasimomentum q belongs to the resonant cou- 
pling line A q = 0, then M q = w q = 1/2. In this case the 
longitudinal optical absorption (|6^) has a split double 
peak at frequencies w+,i q and W-^q, instead of a single 
main peak. The transverse optical absorption (|66|), sim- 
ilarly to the non-resonant case (|66|), has a split double 
peak at frequencies cj +i i q and w_,i q , but its amlitude is 
now of the order of unity. For |q| — > Eq.(|65|) reduces 
to that for an array of noninteracting wires ([23]), and the 
transverse optical conductivity j66| ) vanishes. 

The imaginary part of the ac conductivity a"-, (q, w) 
is calculated within the same approach. Its longitudinal 
component equals 



APPENDIX D. TRIPLE QCB SPECTRUM 

To diagonalize the Hamiltonian (|3^) , we write down 
equations of motion 



[cJ 2 s (qj) - uj 2 ] 0j sq 
Vs^s{qj)^s{qj)— ^0s'(<23V S '(<73)#3s'q = 0, (67) 



Vecj) s {q 3 )uj s (q 3 ) — ^2 $s' )<<V (qj ) 



'js'q 



(68) 



Here j = 1,2, and e is defined by Eq.([L9|). The solutions 
of the set of equations (p7|) - (f38|) have the form: 



.4 



M<li)u s (qj) 



.7 = 1,2,3, 



Substituting this equation into Eqs.(67) and (68), we 
have three equations for constants Af 

A x + A 3 ^F Q3 (lo 2 ) = 0, 
A 2 + A 3 V^F q3 (cu 2 ) = 0, 
A 3 + J2 A j V^F qj (u, 2 ) = 0, 

.7=1,2 



where 



ro 4> 2 s {q)uj 2 (q) 
« „ w s0?) ~ u 2 ' 



(69) 
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Dispersion relations can be obtained from the solvability 
condition for this set of equations 



eF q3 {u 2 ){F qi {uo 2 )+F q2 {u 2 )) = 1. 



(70) 



The function F qs (ui 2 ) has a set of poles at uj 2 — 0J 2 (q), 
s = 1,2,3,... . For uj 2 < uj 2 {q), i.e. within the inter- 
val [0, u>i(q)], F qs (uj 2 ) is positive increasing function. Its 
minimal value F on the interval is reached at u> 2 = and 
does not depend on quasi-momentum q 



F q (0) 



n ^ — ' 



(?) 



(71) 



If the parameter e = rj 2 is smaller than the critical value 



2F 2 '' 



(72) 



then all solutions uj 2 of the characteristic equation are 
positive. When e increases, the lowest QCB mode soft- 
ens and its square frequency vanishes in whole BZ at 
e — e c . For exponential model £(£) = exp(— 1£|), one ob- 
tains e r ~ 1. 



APPENDIX E. TRIPLE RABI OSCILLATIONS 



uj± = ^li v 7 ^ 2 , 
and 0q,± are the corresponding amplitudes. 
Choosing initial conditions 

0i(O)=t0 Ol 9 1 (Q) = cj 9 , 



(73) 



9 2 (0) = 0, 
9 3 (0) = 0, 



9 2 (0) = 0, 
9 3 (0) = 0, 



we obtain for the field amplitudes at the coordinate origin 
0i(O,O;t) = 



sin(w + t) -| sm{LO-t) 

UJ + UJ- 



2(0,0;*) 



j(0,0;t) 



■ sin(wot), 

— sm(cj + r) H sin(w_t) 



-y sm(w *), 

#0 w . , . Wo . , , 

— = — sm(w + £) sin(w_t) 

2v2 L a -'+ W_ 



In the limiting case e <§; 1 these formulas coincide with 
Eqs.(|ill) in subsection III.C. 



The point C(Q/3, 2Q/3) of the BZ is the point of three- 
fold degeneracy, 



qi = 93 = -92 



3' 



wn(Q/3) = wai(2Q/3) - w 3 i(Q/3) = wq. 

Equations of motion at this point in the resonance ap- 
proximation read 



- d 2 
_dt 2 


+ ^0 


9\ - 




h = o, 


- d 2 
_dt 2 


+ ^0 


02 - 




h = o, 


- d 2 

dt 2 


+ ^0 


$3 ~ 


\- y/l4> 2 ul 


(0i + 2 ) = 0, 



where 8j = 0j q . General solution of this system looks as 

9 1 (t)\ ( 1 
2 (t) =0 O -1 e 1 ^* 



7 2l 

?s(*) 








where one of the eigenfrequencies coincides with too , while 
the two others are 
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